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Abstract 

We discuss an application of QCD sum rules to the heavy baryons Aq and £q. The 
predictions for the masses of heavy baryons, residues and Isgur-Wise function are 
presented. The new results on two loop anomalous dimensions of baryonic currents 



00 \ and QCD radiative corrections (two- and three- loop contributions) to the first two 

Wilson coefficients in OPE are explicitly presented. 

m 
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\D '• 1 Introduction 

Baryons with one b quark are very nice systems for application of Heavy Quark Effective 
Theory (HQET)[|I[. It allows one to organize the determination of the properties of 
baryons in an 1/mQ-expansion, the leading term of which gives rise to the spin-flavour 
symmetry. Well known predictions of the HQET are the relations between different hadron 
transition form factors. For example, the six form-factors describing the — > A c electro- 
weak transitions are reduced to one universal Isgur-Wise function £(v ■ v') in the HQ 
limit [0, [3], |4j. However, one still remains with many non-perturbative parameters, such 
as the mass of the baryon, the Isgur-Wise function and the averaged kinetic and magnetic 
energy of the baryon, which should be estimated by using some non-perturbative method 
as sum rules ||, lattice calculation or some potential model. 

In the present paper we review an application of the QCD sum rule method to calculate 



the masses, residues of heavy baryons and baryonic Isgur-Wise functional, 12, 13, 14, 15 



16, 17, IS, IS, ECL 22, 23, 24 



2 Currents, anomalous dimensions and residues 

The currents of the heavy baryon Aq and the heavy quark spin baryon doublet {Sq, £q} 
are associated with the spin-parity quantum numbers j p = + and j p = 1 + for the light 

1 Yakovlev@dipmza.physik.uni-mainz.de 
Invited talk given at the III German- Russian Workshop on Heavy Quark Physics, Dubna, Russia, 
May 20-22, 1996 to appear in the Proceedings 



diquark system with antisymmetric and symmetric flavor structure, respectively. Adding 
the heavy quark to the light quark system, one obtains j p = | + for the Aq baryon and 
the pair of degenerate states j p = | + and j p = | + for the baryons £q and Eg. The 
general structure of the heavy baryon currents has the form (see e.g. [H and refs. therein) 

J = [q iT CVTq']V'Q k e m . (1) 

Here the index T means transposition, C is the charge conjugation matrix with the 
properties C^JC -1 = — 7 M and OyJC -1 = 75, i,j, k are colour indices and r is a matrix 
in flavour space. The effective static field of the heavy quark is denoted by Q. For each of 
the ground state baryon currents there are two independent current components J\ and 

Jai = [q iT CT-f 5 q j ]Q k e ijk , J A2 = [q iT CT l5 -f q j }Q k e ijk , 

■hi = [q iT Criq j ]-^ 5 Q k e ijk , J E2 = [q iT Cr^q j ] ■ ^ 5 Q k e ijk . (2) 

The anomalous dimensions of the heavy baryon currents up to second order was calcu- 
lated recently in [22 and read (for example in the MS scheme with naively anticommuting 

7s) 

7M = -8 (g) + i(16C(2) + 40N F - 796) (^f , (3) 
7A2 = -4 (£) + |(16C(2) + 20N F - 322) [^f , (4) 



7si 



a. 



^(16C(2) + 20iV F - 290) f^V, (5) 
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^ = -KS) + ^ 48 « 2 ) + 8 * + 324) (S) 2 - (6 > 

The anomalous dimensions are ingredient of renormalization group invariant sum rules. 
Numerically effects of 72 are very small and gives about 1 — 2% corrections to the current 
redefinition. 

Let us define residues (i = A, S) of the baryonic currents according to 

(0\J(ri\A Q ) = FaOOt*. (0|J(/*)|£g) = FeOOt*. (7) 

where u is spinor. Then we define bound state energy E R of the baryon R and use the 
formula 

Q 

m b aryon = m Q + E R H hO(l/m Q ), (8) 

m Q 

where mg is the running heavy quark mass at the scale of the same mass and c is a 
coefficient connected with the averaged kinetic and magnetic energy in the baryonic state. 
The coefficient c was recently estimated in [O] by using QCD sum rules. 



3 Correlator of two baryonic currents 

In order to obtain information about the value of the mass of the baryon and its residue 
one considers the correlator of two baryonic currents 

IL(u = k ■ v) = i J d i xe ikx (0\TJ{x)J{0)\0), (9) 



2 



where k^ and 



are the residual momentum and four velocity in = mqv^ + k 



A" 



respectively. P(oS) can be factorized into a spinor dependent piece and scalar function 
P(uj) according to 



n( W ) = r'i±^f'iTr(rf)2Tr(rf)P(^). 
The correlator function Pope(^) satisfies the dispersion relation 



(10) 



Pope (a;) 



p(uj')dio' 



UJ' 



UJ 



+ subtraction, 



where p(uj) = Im(P(uj))/ir is the spectral density. 

Following the standard QCD sum rules method || the correlator is calculated in 
the region — uj ^1 — 2 GeV, including perturbative and non-perturbative effects, where 
non-perturbative effects can be quite important. The non-perturbative effects are taken 
into account with the help of an operator product expansion. The OPE of two point 
correlator was discussed in Refs. [11], |12|) [01 HH- l/ m Q corrections to two and three 



point correlators were calculated in |]19 , [20] , |17| . The leading perterbative term and the 



next-to-leading term in OPE (gluon condensate contribution) give the spectral densities 



L) 



20tt 4 



a s (GG) 

c — OJ. 

32tt 3 



(12) 



Next we consider the radiative corrections to the spectral density of the perturbative con- 
tribution. There are four different graphs contributing to the correlator, which are shown 
in Fig. 1. The fact that all graphs in Fig. 1 have two-point two- loop subgraphs greatly 
simplifies the calculational task. One can first evaluate the respective subgraphs which 
leaves one with a one-loop integration. The first two-loop integration can be performed 
by using algebraic method described in [21 . It is important to note that the results of the 



two-loop integration are polynomials of the external momentum relative to this subgraph. 
Hence, the last integration is really a one-loop one, but the power of one of the propaga- 
tors becomes a non-integer number. The evaluation of results in a Taylor expansion in 
1/e gives |24] 




(13) 



6 | 40C(2) 



61n 2 - 234n 



396 (n-2)s 



p ) \ 45e 2 225e 90 

5(195n 2 - 780n + 1946)((2) - 2200C(3) + 4907n 2 - 18408n + 34352 



2250 



Here we assume that matrix Ti in Eq.(l) is an antisymmetrized product of n Dirac 
matrices: Y = 7^ • ■ ■ j^. This introduces a n- and s-dependence in the QCD corrections 
due to the identities 



hT 



-l) n {D - 2n)r, 7o r 7o 



■l) n sT. 



(14) 
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The correlator is renormalized by the renormalization factor of the baryonic current, which 
is derived in [14], 



P{u) = Z 2 jP rcn (uj) with Zj = 1 + ^^(n 2 - An + 6). (15) 

The multiplication with Z 2 results in the cancelation of the second power in 1/e. The 
first power in 1/e is purely real and hence does not contribute to the spectral density. The 
renormalized spectral density p ven (uj) = Im(P rcn (uj))/7T must in fact be finite and can be 
read immediately off Eq. (]T3|) , 

,5 



1 H r(u/fi) 

7T 



where Po(oo) = and (16) 

20tt 4 



r(u//i) = -(n 2 -4n + 6)ln (— ) + — (60((2) + 38n 2 - 137n + 273). 
3 V 2u ) 45 

The as correction can be seen to depend on the properties of the light-side Dirac matrix 
in the heavy baryon current. The results in the naively anticommuting 75-scheme (AC) 
are 

, , ^ „, ( V \ 2(20C(2) + 91) 
r A1 (^) = 41n(|-) + I y \ (17) 

W16.51 

r^lA - 2., i m + 4(1 ° C(2) + 29) 



2u J 15 



1 12.12 



mWn) = H/n f— ) + 2(6 ° C(2) + 151) . 
v /p; 3 V2cj/ 45 

wii.10 

The coefficient of the logarithmic term ln(2a;//x) in Eq. ([IB] ) coincides with twice the one- 
loop anomalous dimension given in Eq. (|i~5[ ), as expected. The results for the Baryons Ai 
and A2 in the 't Hooft-Veltman 75-scheme (HV) differ from those presented above. But 
we know that currents in different schemes should be connected by finite renormalization 
factors Z, Jac = ZJhv- Such factors was recently derived from two-loop anomalous 
dimensions of baryonic currents f22j. They read 



Z A1 = 1-^ and Za2 = 1-^. (18) 
Sir oil 

By multiplying the results in the 't Hooft-Veltman scheme by Z 2 we obtain the same 

results as in the naively anticommuting 75-scheme. 

The results show that the as-corrections amount to about 100%, which makes per- 

turbative QCD radiative "corrections" very important in QCD sum rules. We used 

fj, — 2uj w 1 GeV and a>s(p) = 0.3. 



4 Contribution of quark condensate 

Next, we consider the contribution of the quark condensate, which appears in the nondi- 
agonal sum rules. The leading and next-to-leading spectral density are 

, ( „H-<^ * M = 2 (i_£)M=3. (19) 
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The radiative corrections to the quark condensate term can be quite important, and we 
take it into account. There are 8 different graphs contributing to the correlator, which 
are shown in Fig. 2. We used the same method as we used for the perturbative term. The 



correlator must be renormalized by the renormalization factor of the baryonic current [14] 
and quark condensate, 

P(u) = Z A Zj 2 Z qq P™{oo) with Z, h = l + ^^(n 2 -4n + 6), (20) 

47T6 

Zj 2 = 1 + ^=f((n - 2 + s) 2 + 2), Z- qq = 1 + |j^(-3). (21) 

This procedure cancels the second power in 1/e while the first power in 1/e is purely real 
and does not contribute to the spectral density. For the renormalized spectral density we 
get 

= _«^[ 1 + ^(* (an ._8 B + 7 + 2(B _ 2)4)bl (|. 

+ -(8n 2 - 28n + 37 + 8ns - Us + 8C(2)))] . (22) 

Then we proceed with the usual QCD sum rules analysis and equate the theoretical result 
for Pope(^) with the dispersion integral over the hadron states saturated by the lowest 
lying state with the bound state energy Er plus exited states and continuum. So the 
phenomenological part of the sum rule is given by the spectral density p p h(^) = p rC s(^) + 
Pcont(^), where the contribution p res of the low-lying baryon state is p r es(^) = ^-5(tu — 
Er). As is usual for the contribution of excited states and continuum contributions we 
assume hadron-parton duality and take p CO nt(^) = 0(u — Ec)p(uj), where p is the result 
of the OPE calculations. We apply the Borel transformation 



B T = lim — — n, —lu — » oo, T = —lu/ti fixed (23) 

T(n) \ doj ) 

and finally obtain the sum rule 

1 rEr 

L F i e -E R /T = /• p {jy-u'IT d j + BP pc (T). (24) 

2 J o 



5 Numerical results 

Let us discuss the sum rule analysis. First, we analyse the dependence of the bound state 
energy Er as a function of the energy of continuum Eq and the Borel parameter T in a 
wide region of these arguments and try to find regions of stability in T and E c . Second, we 
fix the energy of continuum E c = E^! est by requiring that Er(T) should have the highest 
possible stability in its dependence on the Borel parameter T. We take into account the 
as-correction to the spectral density and obtain 

E A = 0.9 ± 0.2 GeV and E^ est = 1.2 ± 0.2 GeV. (25) 

The analysis of sum rules at E\ = 0.9 GeV gives for the A-residue 

F A = (0.03 - 0.04) GeV 3 and E c = 1.0 - 1.4 GeV. (26) 
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Exp. |26| and 

Pot monpl \?T\ 


[141 


1171 
1 1 ' \ 




L 


N.L.O. 


E C (A) 




1.20 
1.46 


1.2 ±0.1 
1.4 ±0.1 


1.20 ±0.15 
1.30 ±0.15 


1.2 ±0.2 
1.4 ±0.2 


1.2 ±0.2 
1.4 ±0.2 


Er[A) 
Er(X) 
Diff 

LJ ILL. 


U.o4U(^l.U4UJ 

1.050(1.250) 
n i qo 


U. m 
0.99 


n n -L n 1 

u.y ± u.i 


U. ( y ± U.UO 

0.96 ±0.05 
n 1 7 

U.I / 


U.o ± U.Z 
1.02 ±0.2 

n 9 


n no i f\ o 

u.y/ ± U.Z 
1.12 ±0.2 
n 9 


A(A) 
ACE) 


0.460 
0.405 








0.40 ±0.05 
38 ± 05 


0.32 ±0.05 
32 ± 05 


Fa 




2.3 ±0.5 
3.5 ±0.6 


2.5 ±0.5 
4.0 ±0.5 


1.7 ±0.6 
4.1 ±0.6 


2.4 ±0.5 
3.4 ±0.5 


3.2 ±0.5 
5.2 ±0.5 



Table 1: All energies like Ec, Er and A are given in GeV, the residues in 10" 2 GeV 3 . The 
value of the Borel parameter is T = 0.6 GeV. We use m& = 4800 MeV (resp. 4600 MeV) 
in the first colomn. 



Doing the same analysis for the E baryon and taking into account the a^-corrections we 
obtain 

£ s = 1.1 ± 0.2 GeV and E^ est = 1.4 ± 0.2 GeV. (27) 



F s = (0.045 - 0.055) GeV 6 and E b c cst = 1.4 ± 0.2 GeV. 



:28) 



All the results are summarized in the Table 1 , where we compare our results with the 
leading order results obtained in the |14, K% |19[ and some experimental values. 



6 The Isgur-Wise function 

Next, we consider semileptonic transition A& — > A c . The matrix elements of the weak 
current for A^ — > A c at the leading order 1 / rriQ are determined only by Isgur-Wise function 

< AelcT&IA, >= t(y)u c Tu b . (29) 

where y = v j ■ v c and T is some gamma matrix. The Isgur-Wise function was calculated 
in Ref.fl5| by using three-point QCD sum rules. The slope of Isgur-Wise function at y=l 



is 

p 2 = -1.15 ±0.2. (30) 

The uncertainty is connected mainly with an assumption about continuum model. We 
have found that the shape of the Isgur-Wise function nearly coincides with the ansatz 
formula 

m = -^-^pW - 1)^], (si) 



with the slope in Eq. fl3"0|) . Fig. 3 shows the Isgur-Wise function at different energies of 
the continuum (E c = 1.0, 1.2, 1. AGeV). 1/toq corections to A& — > A c transition were 
discussed recently in [pD| . The important 1/mg effects to the decay come only from the 
weak current expansion. 
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7 Conclusions 



We have calculated one and two-loop anomalous dimensions of baryonic currents. We have 
studied the expansion of the correlator of two heavy baryon currents at small Euclidian 
distances. The radiative corrections to the first two Wilson coefficients are calculated. 
The heavy baryons sum rules in order are derived. The predictions for the mass, 
residues , the slope and shape of Isgur-Wise function are presented. 
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(0) 




(3) (4) 

Figure 1 



Figure 1: Two-loop and three-loop contributions to the correlator of two heavy baryon 
currents. (0) two loop lowest order contribution, (l)-(4) three loop 0(cts) contributions. 
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(6) (7) 



Figure 1 



(2) 




(8) 



Figure 2: One- loop and two-loop contributions to the nondiagonal correlator of two heavy 
baryon currents. (0) one loop lowest order contribution, (l)-(8) two loop 0(as) contri- 
butions. 
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Figure 3: The Isgur-Wise function. 
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